Physics counterpart of the VT non-hermitian tight-binding chain 
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We explore an alternative way of finding the link between a VT non-Hermitian Hamiltonian and 
a Hermitian one. Based on the analysis of the scattering problem for an imaginary potential and its 
time reversal process, it is shown that any real-energy eigenstate of a VT tight-binding lattice with 
on-site imaginary potentials shares the same wave function with a resonant transmission state of 
the corresponding Hermitian lattice embedded in a chain. It indicates that the VT eigenstate of a 
VT non-Hermitian Hamiltonian has connection to the resonance transmission state of the extended 
Hermitian Hamiltonian. 
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I. INTRODUCTION 

Imaginary potential usually appears in a system to de- 
scribe physical processes phenomenologically due to its 
simplicity, which has been investigated under the non- 
Hermitian quantum mechanics framework [U-Hoj . To dis- 
cuss and explore the usefulness of the imaginary poten- 
tial, one has to be able to establish a correspondence 
between a non-Hermitian system and a real physical sys- 
tem in an analytically exact manner. The effort to es- 
tablish a parity-time {VT) symmetric quantum theory 
as a complex extension of the conventional quantum me- 
chanics jll>tl7j was stimulated by the discovery that a 
non-Hermitian Hamiltonian with simultaneous VT sym- 
metry has an entirely real quantum-mechanical energy 
spectrum [l8| and has profound theoretical and method- 
ological implications. 

When speaking of physical significance of a non- 
Hermitian Hamiltonian, it is implicitly assumed that 
there exists another Hermitian Hamiltonian which shares 
the complete or partial spectrum with it. Thus one of 
the ways of extracting the physical meaning of a pseudo- 
Hermitian Hamiltonian having a rea l sp ectrum is to 
seek for its Hermitian counterparts [a-Hoj. The metric- 
operator theory outlined in [ll| provides a mapping of 
such a pseudo-Hermitian Hamiltonian to an equivalent 
Hermitian Hamiltonian. However, the obtained equiva- 
lent Hermitian Hamiltonian is usually quite complicated 
(III HH and cannot be judged whether it describes real 
physics or is just unrealistic mathematical object. 

In this paper, we try to find an alternative way to 
establish the connection between a pseudo-Hermitian 
Hamiltonian and a physics system. We consider a simple 
class of discrete systems, which are originally exploited 
to describe the solid-state system in condensed matter 
physics. In such systems, the imaginary potential usu- 
ally appears as source or sink, acting as connection to the 
outer world. In this sense, the eigenstates of an unbroken 
VT non-Hermitian Hamiltonian seem to be the dynami- 



cal equilibrium states of an open system. The strategy of 
this paper is to seek the ways of analytical continuation 
of the eigenfunctions of a VT non-Hermitian Hamilto- 
nian into the stationary scattering states of an extended 
Hermitian system. It indicates that the VT eigenstate of 
a VT non-Hermitian Hamiltonian has connection to the 
resonance transmission state of the extended Hermitian 
Hamiltonian. 

This paper is organized as follows, in Sec. II, Scat- 
tering problem of imaginary potential. In Sec. Ill, the 
connection between a resonant transmission state and a 
VT symmetry eigenstate is established. In Sec. IV, the 
illustrative examples are presented to demonstrate the 
main idea of this paper. Sec. V is the summary and 
discussion. 



II. SCATTERING PROBLEM OF IMAGINARY 
POTENTIAL 

The first problem we address in our search for physi- 
cally meaningful VT non-Hermitian Hamiltonian is how 
to associate the individual imaginary potential with the 
Hermitian sub-network. Recently the formal theory of 
scattering for complex potentials in one dimension con- 
tinuous system has been constructed (for review see Ref. 
20] and references therein). To establish such a cor- 
respondence in a discrete system, we consider a simple 
model described by a non-Hermitian Hamiltonian H 1 . 
It is a tight-binding chain with uniform nearest neigh- 
bor hopping integrals and an additional imaginary on- 
site potential on a site of a semi-infinite chain, which can 
be written as follows: 



H 1 — Hiy + H g + -ff su b, 
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FIG. 1: (Color online) Schematic illustration of configurations 
of the typical tight-binging networks with imaginary poten- 
tial and their Hermitian counterparts, (a) the scattering pro- 
cess of an absorptive (absorbing) imaginary on-site potential, 
where H su h is an arbitrary Hermitian sub-network, (b) the 
corresponding Hermitian counterpart network of (a), which 
ensures the same wave function as (a) within the common 
region. The system consists of system (a) and an attaching 
lead, (c) and (d) describe the situations of (a) and (b) un- 
der the time reversal operations. Note that (c) represents the 
scattering process of a new network different from (a), while 
(d) corresponds to the time reversal scattering process of (c) 
in the same system. 

where a[ is the creation operator of the boson (or 
fermion) at the /th site, the tunneling strengths and 
imaginary potential are denoted by J, <?, and —ij (7 > 0). 
We separate the Hamiltonian into severval parts to char- 
acterize the configuration of the system. A sketch of such 
a system is given in Fig. 1(a). Here Hz 7 is a semi- infinite 
uniform chain with one potential —ij at the edge, H g 
represents the coupling between this chain and an arbi- 
trary sub-network described by a Hermitian Hamiltonian 
i/ su b- In this sense, the conclusion will be obtained is ap- 
plicable for a large class of systems. 

To investigate the role of the imaginary potential in a 
discrete system, we will be concerned with the scattering 
problem of such system: an incident plane wave e 1 ^ or a 
broad wave packet comes from leftmost and is reflected 
and transmitted at the imaginary potential. The pro- 
cess can be represented by the wave function f\ (j) a J |0) 
(j G (-00, N s }) with 

/i(i<0) = e i ^+f 1 e-*, (2) 

where r\ represents the reflection amplitude. The explicit 
form of wave function f\ {j > 0) within the sub-network 
depends on the the structure of -ff S ub- General speaking, 
the solution of r\ can not be obtained exactly even the 
explicit form of -ff su b is given. However, we will see that 
it does not affect the conclusion below. 

In the basis < |0) \j € (—00, N s ] [, The Schrodinger 



equation has the explicit form 

-J/i(j-l)- Jfi(j + l)=Efi(j), (j<0) (3) 

-J/i(-i)-<?/i(i) = (£ + n)/i(o), 

-fl/iM+E" 1 ! «a/i (*) = -B/i (1) 
£^>«/i» = £/i(j), (ie [2.JVJ) 

within all the regions. We will show that, such a scatter- 
ing process can occur in a Hermitian system. 

It has been well known that an imaginary potential, 
by means of an effective interaction, can serve as a re- 
duced description for the outer world of an open system. 
Along this line, we consider a similar lattice system to 
Hj, described by a Hermitian Hamiltonian Hy. In this 
network, the imaginary potential is replaced by a real po- 
tential V, and a semi-infinite chain is added, which acts 
as the complementary subspace or outer word if H 7 is 
regarded as an open system. The corresponding Hamil- 
tonian has the form 

Hy — Hiv + H g + H su b + H v + Hi v (4) 

-1 

Hiv = -J 51 (ajai+i + H.c.) + Vajao 

/ — — OO 

H v = -v {a^b v + b\ao^ 

00 

H lv = -jJ2(b]bi +1 +R.c^ . 

1=1 

where b) is also the creation operator of the boson (or 
fermion) at the Ith site. Hiv and Hi 7 are two semi- 
infinite uniform chains with real potential V at the joint 
point. 

A sketch of such a system is given in Fig. 1(b). 
Note that the two systems have assigned the same sub- 
network, j € (—oo,0], which is referred as the com- 
mon region. The corresponding scattering wave func- 
tions within the two semi-infinite chains are f 2 {j) cij |0) 
and f(j)bt |0) with 

h (j < 0) = e tt ' + r 2 e-^ (5) 
f(j>0) = te k K 

In the basis {aj |0) \j € (-00, N s ],b] |0) \j € (1, +00]}, 
the Schrodinger equations are 

-Jf2(j-l)-Jf2(j + l) = Ef 2 (j), (j<0), (6) 
-J/ 2 (-1) - gf 2 (1) - uf(l) = (E-V) h (0) , 

N„ 

-ah (0) + 2 (0 = -E/2(l), 

1=1 

Y,Kijh(i) = Ef 2 (j), (je[2,N s }), 

i=l 
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and 

-vf 2 (0)-Jf(2) = Ef(l), (7) 
-Jjy-l)- Jf(j + 1)= Ef(j), (j>l). 

We can see that the Eqs. <j3j> and (|6|) within the common 
region have the same form except the potentials at the 
Oth site. For the the same incident plane wave e* fcj ', we 
have 

E = -2Jcosk (8) 

and furthermore, one can find that under the conditions 

v 2 sin k = 7 J, v 2 cos k = VJ, (9) 

the solutions for r\ and ri are identical. The above 
equivalent conditions can also be given in the energy- 
dependent form 
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FIG. 2: (Color online) Schematic illustration of configurations 
of a non-Hermitian tight-binging network with VT symmetry 
and its Hermitian counterpart. They are constructed based 
on the building blocks represented in Fig. 1. It indicates 
that an eigenstate (stationary state) of (a) corresponds to a 
resonant transmission state of (b). 
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where f2 = V '4 J 2 — E 2 . Then the wave functions @ and 
([5]) within j £ [— oo, N s ] are the same. This indicates 
that for the incident plane wave e 1 ^ , the imaginary po- 
tential can be regarded as a semi-infinite chain for wave 
escaping. It is worth to point out that it is a conditional 
equivalence, which is only for the specific state. This 
equivalence is the building block for the investigation of 
this paper. 

Now we consider the relevant situations derived from 
the obtained results. We are interested in the case when 
the imaginary potential is source like. Actually, apply- 
ing time-reversal operation on the above scattering pro- 
cesses, i.e., taking complex conjugation for the Eqs. ([T]) 
and ([3]), the corresponding time reversal solutions can 
be obtained, which are illustrated in Fig. l(c, d). The 
Fig. 1(c) shows that the time reversal solution is for the 
new system = H 7 (7 — > —7). Nevertheless, its coun- 
terpart Hy is invariant under time reversal. The Fig. 
1(d) illustrates the corresponding time reversal process 
of that in Fig. 1(b). Based on the processes in real phys- 
ical systems illustrated in Fig. l(b, c), the physics of the 
imaginary potential becomes clear: —ij acts as a drain 
lead, while 17 acts as a source lead associated with an 
incoming plane wave. Although this is not a surprising 
result, we still verify it explicitly in a strict manner and 
will apply it to a VT non-Hermitian system. 



III. RESONANT TRANSMISSION CONDITION 
AND VT SYMMETRY 

The above result is essentially about the stationary 
state for the infinite non-Hcrmitian system. Intuitively, 
a source and drain could produce a stationary state in a 



finite system between them when the gain rate equals 
to the loss rate, or in an open system with injecting 
sources and absorbing sinks. The conceptual framework 
is required to substantiate this idea. In this section, we 
consider the stationary state of a non-Hermitian system 
based on the obtained scattering solutions of the imagi- 
nary potential. 

Although we cannot get the explicit solution about the 
reflection amplitude r, the corresponding time-reversal 
process illustrated in Fig. 1(c) exhibits the facts: for an 
incident plane wave with amplitude 1, the reflected am- 
plitude from the absorptive potential — 27 is r, while for 
an incident plane wave with amplitude r*, the reflected 
amplitude from the source potential 17 is just 1. The 
fact \r\ — \r*\ indicates that if we combine the build- 
ing blocks Fig. 1(a) and (c) to construct a finite network 
with the geometry illustrated in Fig. 2(a), the stationary 
state may be formed in the manner: a wave coming from 
the sources and send back its time-reversed version. On 
the other hand, such a configuration has the VT sym- 
metry spontaneously, which has been shown to process 
real-energy eigenstate under certain condition. In the 
aid of its Hermitian counterpart shown in Fig. 2(b), one 
can find that in the case of the stationary state being 
formed, it just corresponds to the resonant transmission. 
Therefore the existence of the real-energy eigenstate of a 
VT symmetric system is connected to the occurrence of 
resonant transmission in a Hermitian system. It follows 
that we can find an alternative Hermitian counterpart 
to a VT symmetric Hamiltonian in the sense that they 
share the same eigenfunction within the common region. 
This should be the way more directly associated with the 
physics of the VT symmetric system. In the following 
sections, we will study the formation of the VT symmet- 
rical state by dealing with the more tractable models. 
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Note the region {aj |0) , le (l,iV + 2iV s )} is regarded 
as the common region of the two models. In the follow- 
ing, we present the analytical results in the framework of 
Bethe ansatz for the two models in order to perform a 
comprehensive study. 
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FIG. 3: (Color online) Schematic illustrations of tight-binging 
network with VT symmetry and its Hermitian counterpart, 
where -ff su b is a simple chain of length N s . Exact Bethe ansatz 
solution shows that the eigenstate of (a) accords to the equal- 
energy resonant scattering state of (b) under the condition 
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IV. ILLUSTRATIVE EXAMPLES 

In this section, we investigate a simple exactly solvable 
system to illustrate the main idea of this paper. In or- 
der to exemplify the above mentioned analysis of relating 
the stationary states of a non-Hermitian VT symmetric 
Hamiltonian and a Hermitian one, we take H su ^ to be 
the simplest network: a uniform chain. Then the sample 
Hamiltonian has the form 



N+2N S -1 



(ii) 



+il a N s +i a N s +i - iia' ! N+Ns a N+Ns 



which is sketched in Fig. 3(a). It has VT symmetry, 
i.e., H^T — VTHy^VT — H^, which has been studied 
systematically in the case with zero N e . Here V and 
T represent the space-reflection operator (or parity op- 
erator) and the time-reversal operator respectively. The 
corresponding Hermitian Hamiltonian reads 



N+2N S -1 

J a}ai +1 - jJ2 b l b i±i (I 2 ) 

i=i ;=±i 

-v Uijv 8+ i&-i + ajv s +jv & i) + Hc 

+y (< 1 N s +i a N 3 +i + ajvs+AraJVs+A') , 

which is sketched in Fig. 3(b). It is a V symmetric 
system, i.e., [P,H] = 0, which has been studied in the 
framework of Bethe ansatz for case with zero V [2l| . The 
effects of V and Ton a discrete system are 



TiT 



-i, Va\V 



a N+2N s +2-n ^{V 



(13) 



A. VT chain 

According to the ■pT-symmetric quantum mechanics 
[l6j |. system H 1 ^ l can be further classified to be either 
unbroken VT symmetry or broken VT symmetry, which 
depends on the symmetry of the eigenstates in dif- 
ferent regions of 7. The time-independent Schrodinger 
equation is 



(14) 

with corresponding eigenvalue e7. The system is unbro- 
ken VT symmetry if all the eigenfunctions have VT sym- 
metry 



VT^l) 



(15) 



and all the corresponding eigenvalues are real simulta- 
neously. This classification depends on the value of the 
parameter 7. Beyond the unbroken VT symmetric re- 
gion the system is broken VT symmetry, where Eq. (|15[) 
does not hold for all the eigenfunctions and the eigenval- 
ues of broken VT symmetric eigenfunctions are complex. 
We denote the single-particle eigenfunction in the form 

W> = £# (0*J|o>- 

According to Bethe ansatz method, the eigenstates can 
be in the form of 



C L e ik i + D L e~ tk J, j e [1, N s + 1] 
fl = { Ae ik * + Be~ lk \ [N. + 1, N s + N] . (16) 
C R e ikj + /)/,■< ' ? '-' , [N S +N,N + 2N S ] 

The coefficients { A, B, Cl(.r) > -Cl(-R)} an£ l the quasi mo- 
mentum k are to be determined by matching conditions 



AH? + o + )=A 7 (j + o-) 

(j = N s + l,N s + N) 



(17) 



and the corresponding Schrodinger equations of j (j ^ 
N s + 1, N + N s ) in the center of the system, 



(18) 



the Schrodinger equations of j for the sites j = N s + 1 
and N, + N, 



-jP k (N s +2)-Jf k {N s ) 

(4-n)/ fe 7 (iV s + i), 

-Jfl (N s + N + 1) - Jf k (N s +N- 1) 



(19) 
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and for the edges of the system j = 1 and N + 2N S 

-jP k {2)=elP k {l) 
-J ft (N + 2N S - 1) = elfl (N + 2N S ) . 



These lead to the equation of k 



2 2 

"7 Xk 



Jk(N-l) 



ik(N+2N B + l) 



3 — tfc(iV-l) 



-ik(N+2N s + l) 



(20) 



(21) 



where 



Jk(N s + l) 



-ik(N s + l) 



Xk 



-ik 



(22) 



The solutions of k can be classified in two categories: 
physical and unphysical states, which correspond to real 
and complex k, respectively. The corresponding energy 
are real or complex and in the form 



-J(e 



ik 



-ik 



(23) 



A straightforward algebra shows that there are at least 
(N — 1) solutions of real k for arbitrary 7/ J. In this 
paper, we only focus on the physics counterparts of these 
states rather than the detailed form of the solutions. 



B. Hermitian counterpart Hvv 



For the Hamiltonian Eq. (|12l) . it has V and T sym- 
metry simultaneously, i.e., VHyv'P — THvvT = Hvv- 
Nevertheless, for a scattering state, a plane wave comes 
from the leftmost, the V and T symmetry are broken. 
We will show that under certain conditions, the wave 
function within the common region of and Hvv has 
VT symmetry. 

We can set the scattering wave function in the form of 

C| e ifei + D s L e- ikj , [l,iV s + l] 
fZ = { A s e ik i + B s e-^, [N S + 1,N S + N] , 
C R e ik i + D s R re- ikj , [N s + N, N + 2N S 



and 



fk 



e tkj _|_ re 



-ikj j £ (-00, -1] 
[l,+cx>) 



(24) 



(25) 



where f k represents the one within the common region, 
while f% represents the one in the two leads. The re- 
flection and transmission amplitudes r, t, coefficients 

(a s , B s , C^ R y ^1,(^)1 and the quasi momentum k are 
to be determined by matching conditions 



(j = N s + l,N + N s ) 
/7 (N s + 1) = f k (N s + 1) 
]l (N s + N) — ft (N s + N) 



(26) 



and the corresponding Schrodinger equations of j in the 
center of the system, 



J/ni + l)- J/ni-l) = elfXV) (27) 

-7 ■ 



-J/T(i+l)-J/7(i-l) = elfj(j) 



Schrodinger equations of j for the connection sites N s + 
1,JV, + JV, 

_ v /7(_l)-J/^(JV. + 2) (28) 
-JfV(N s ) = (el-V)flf(N s + l), 
-vfJ(l)-JfV {Ns + N + 1) 

-JfZ (N s +N -l) = (el-V) & (N, + N) , 
and for the edges of the system j — I. N + 2N S , 

~JfZ(2) - elfZ(l) (29) 
-JfZ(N + 2N s -l) = elfZ(N + 2N s ). 

The solution for reflection amplitude r is given, after a 
straightforward algebra, by 



2w 2 £sin (k) sin 2 [k (N — 1)] 
J 2 sin 2 (jfe) - J 2 e sin 2 [k (N - 1)] 



(30) 



where 



r= y 1 sm[fciV] ^ 2 lk sm[kN s ] 

5 J sin [fc (N — 1)] J 2 sin [fe (iV s + 1)] ' 1 j 

We are interested in the resonant transmission state, 
which should be relevant to the VT eigenstate of H^, 
according to the above time reversal analysis. For r = 
we have 



2w 2 £sin(fc)sin 2 [k (N — 1)] 
= J 2 sin 2 (k) - J 2 £ 2 sin 2 [jfe (N - 1)] 



(32) 



The analysis in Sec. II indicates that, under the con- 
ditions (|9]) , there should be a resonant transmission state 
corresponding to the eigenstate of (|11|) . In fact, substi- 
tuting ([9]) into (|32|) . one can obtain the 

7 2 sin 2 [k (N s + 1)] sin [k (1 - N)] (33) 
= J 2 sin 2 (jfe) sin [jfe (N + 2N S + 1)] , 

which is just the reduced form of the Eq. (|2Tj) for the real 
quasi momentum k. It can be shown exactly that there 
are at least N — 1 real solutions for the Eq. (|3"3")1 . This 
exhibits the connection between the two models (fTTj) and 
(fT2|) predicted by the above mentioned analysis. Accord- 
ingly, we have 



C. 



L(R) 



D 



L(R) 



C 



L(R) 



D 



L(R) 



A_ 



B 
~B~ S 



(34) 
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which shows that both functions within the common 
region are identical, i.e., the scattering wave function 

(r,t, A S ,B S , C^ R y is the analytical continua- 

tion of the one {A, B, C L ( R) , D L ( R) ). 

Thus two eigenstates of the Hamiltonian (jlll) belong 
to the resonant transmission scattering states of two 
different systems (with different V and v). In this 
sense, nonorthogonality of the eigenstates of a Pseudo- 
Hermitian Hamiltonian is obvious. In general the norm 
of a wave function is conserved in a Hcrmitian system, 
while the norm of many components of the wave function 
is not. These wave components may be associated with 
a truncated basis set, or a subspace of the full Hilbert 
space. The common region of and Hvv is a con- 
crete example in the discrete system to demonstrate this 
fact. 



V. CONCLUSION AND DISCUSSION 

In conclusion, we have presented an alternative way 
of finding the link between a VT non- Hcrmitian Hamil- 
tonian and a Hermitian one, based on the analysis of 
the scattering problem for an imaginary potential. We 
have found that the on-site absorptive imaginary poten- 



tial can be equivalent to an attached semi-infinite chain 
as a drain with respect to a specific wave scattering prob- 
lem. Applying this result and its extension to the time 
reversal process on the VT non-Hcrmitian system, the 
eigenstate problem is connected to that of resonant trans- 
mission problem in the corresponding Hcrmitian system. 
It is shown that any real-energy eigenstate of a VT tight- 
binding lattice with on-site imaginary potentials shares 
the same wave function with a resonant transmission 
state of the corresponding Hermitian lattice embedded 
in a chain. 

It is not surprising that the VT eigenstate has connec- 
tion to the resonance transmission state of the extended 
system. In general, for an infinite system with parity 
(V) symmetry, it also has the time (T) reversal symme- 
try. An arbitrary scattering state, as an eigenstate of 
the system, probably breaks the V, T and also VT sym- 
metries simultaneously. Interestingly, in the case of the 
resonance transmission, the corresponding wave function 
possesses VT symmetry. Such wave functions can be the 
eigenstates of the VT non-Hermitian system with the 
real eigenvalues. 

We acknowledge the support of the CNSF (Grants No. 
10874091 and No. 2006CB921205). 
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